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for  Second  Order  Elliptic  Problems"^ 
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The  research  results  reported  in  the  present  paper  represent  a 
substantial  advance  in  the  direction  of  providing  more  efficient,  cost¬ 
saving  techniques  for  solving  a  wide  class  of  commonly  occurring  two- 
dimensional  boundary  value  problems.  In  previous  papers  ([5],  [6]),  it 
has  been  shown  that  it  is  possible  to  dramatically  reduce  the  cost  of 
solving  two-dimensional  problems  by  amalgating  three  formerly  disparate 
problem-solving  tools,  namely: 

1.  Computer  graphics  (visual  feedback) 

2.  Numerical  analysis  (scientific  computing  software) 

3.  Qualitative  information  (the  analyst's  experience  and 
insight,  and  "weak"  mathematical  theorems). 

More  specifically,  in  [5],  Gordon  and  Hall  pointed  out  (via  examples) 
the  practical  utility  of  such  an  amalgamation.  The  problems  considered 
therein  were,  however,  restricted  to  elliptic  boundary  value  problems  sub¬ 
ject  to  Dirlchlet  boundary  conditions,  i.e.,  problems  in  which  the  function 
values  are  specified  on  the  perimeter  of  the  domain.  That  paper,  as  well 
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as  [6],  focused  on  the  Issue  of  contrasting  the  ustial  way  of  Initializing 
an  Iterative  numerical  solution  method  with  the  proposed  new  technique 
which  uses  the  so-called  "blending-function  methods"  of  Interpolation  to 
a  exactly  match  the  boundary  conditions. 

As  one  would  intuitively  expect,  starting  with  what  literally  "looks 
like"  (computer  graphics)  a  good  approximation  reduces  the  computation 
(numerical  analysis)  time  very  substantially.  If,  In  addition,  an  analyst 
Is  provided  a  mechanism  for  quantifying  his  experience-based  knowledge 
(qualitative  information)  of  the  particular  class  of  problems  under  study, 
the  "exact"  solution  is  almost  In  hand. 

The  Gordon/Kelly  paper  [6]  extends  these  early  results,  involving  only 
Dlrlchlet  boundary  conditions,  to  the  rather  general  problem  of  satisfying 

"mixed  linear  boundary  conditions,"  l.e.,  boundary  conditions  of  the  form: 

3F 

oF  +  0  ■  g.  The  boundary  conditions  are,  however,  assumed  to  be  "consis¬ 

tent."  By  this  is  meant  that,  at  the  comers  of  the  region,  the  boundary 
conditions  from  either  side  "match." 

The  attached  paper  addresses  the  problem  of  inconsistently  specified 
boundary  conditions.  In  the  simplest  Instance  of  Dlrlchlet  conditions,  this 
means  that  the  function  values  do  not  match  at  the  comers.  Herein,  we  show 
how  to  actually  construct  bivariate  functions  which  exactly  match  the  above 
type-’inlxed  linear  boundary  conditions,  even  when  they  are  Inconsistently 
specified,  cf.  Section  III.  Moreover,  software  has  been  developed  which 
numerically  performs  the  necessary  operator  multiplications  and  constructs 
the  singular  functions  needed  to  accommodate  such  Inconsistent  boundary 
conditions.  This  software  will  soon  be  available  for  general  distribution.^ 


Contract  Number  N00014-80-C-07H 
Identification  Number  NR  044-56' 


EXACT  A  PRIORI  HATCHING  OF  HIXFD  BOl/NDARY  CONDITIONS 


This  paper  Is  to  appear  In  the 
Prooeedinga  of  the  Foupth  Inter- 
S’j^oaiwi  on  Cerrputer 
Methoda  for  Partial  Differential 
Bguationa  to  be  held  on  June  30  - 
July  2,  1981,  at  Uhlgh  University. 


FOR  SECOND  ORDER  ELLIPTIC  PROBLEMS 

Wtlllaa  J.  Gordon 
Linda  C.  Thiel 

Hatheaatlcal  Sciences  Departawnt 
Drexel  University 
Philadelphia,  PA  19104 


In  this  paper  we  consider  Che  problen  of  con¬ 
structing  bivariate  functions  which  exactly  match 
boundary  conditions  of  the  general  form  oF  B(3F/3n) 
-  «(».y)  on  the  perimeter  of  the  unit  square.  The 
reason  for  wishing  to  do  this  Is  that  substantial 
savings  In  computation  time  can  be  realized  In  the 
subsequent  solution  of  the  discretized  boundary  value 
problem:  If  an  Iterative  method  Is  used  to  solve  the 
discretized  problem,  beginning  with  a  good  Initial 
approximation  can  dramatically  reduce  the  number  of 
Iterations  required  to  achieve  convergence;  If  a  di¬ 
rect  solution  method  is  used,  a  specified  accuracy 
can  be  achieved  with  far  fewer  algebraic  unknowns. 
Although  attention  herein  Is  restricted  to  rectangu¬ 
lar  regions,  the  techniques  developed  can  be 
Straightforwardly  extended  to  any  rectangular 
polygon.  The  Interpolation  techniques  which  we  de¬ 
velop  for  axact  boundary  matching  are  Illustrated  by 
several  exa.iples  which  are  accompanied  by  perspective 
views  of  their  graphs,  and  by  contour  plots. 

I.  Background  and  Introduction 

To  set  the  stage  for  the  main  Ideas  of  this 
paper,  wc  begin  by  displaying  the  familiar  bllinearly 
blended  Interpplant  to  Diri thlet  boundiru  aenditiona 
on  the  perimeter  of  the  unit  square  S  •  [0,l]x(0,l]. 
Let  F(x,y)  be  a  supposed  primitive  function  from 
which  the  boundary  conditions  are  extracted.  Then, 
the  aynthetio  funation  l'(x,y)  which  »re  construct  via 
transflnlte  ("blending  function")  Interpolation  Is 
given  by  the  fallowing  (cf.,  [1],  (2],  [3]): 

0(x,y)  -  (l-x)K(O.y)  +  xF(l.y) 

(l-y)F(x,0)  +  yF(x,l)  ,  , 

-  (l-x)<l-y)F(0.0)  -  (l-x)yF(0,l) 

-  s(l-y)F(l,0)  -  xyF(l,l). 

It  la  easy  to  verify  that  U  •  F  on  3S. 

As  an  example  of  the  use  of  (1.1),  consider  e 
primitive  function  F  whose  values  along  the  four 
edges  of  the  unit  square  are: 

F(0,y)  -  sln(2sy)  ♦  1.5  F(l.y)  •  y^(y-l)  +  .5 


F(x,0)  -  1.5 


F(x,l)  -  (x-1)^  +  .5. 


The  function  U  which  Interpolates  these  boundary  con¬ 
ditions  is  given  by  (1.1)  as: 

0(s,y)  -  (l-x)sln(2i»y)  +  xy(y*-y-2) 

j  (1.3 

♦  s  (2y-l)  +  1.5. 

The  graph  of  U  Is  shown  In  Fig.  la,  and  s  contour 
plot  of  U  Is  given  tn  Fig.  lb. 


This  work  was  supported  hv  the  U.S.  Office  of  Naval 
Research  and  the  Air  force  office  of  Scientific 
Research  under  contract  N00i)l 4-80-0-07 16  to  Drexel 
University. 


Bivariate  and  higher  dlaienslonal  Interpolation  Is 
most  easily  discussed  In  the  formalism  of  projection 
operators  ([1],  [2],  [5],  [6]).  For  Instance,  the 
above  expression  (1.1)  for  U  can  be  more  succinctly 
expressed  as  the  Boolean  am  of  the  two  elementary 
projectors  and  P^  given  by: 

f,lF]  -  (l-x)F(0,y)  +  xF(l,y) 

*  (1.4) 

Py[Fl  -  (l-y)F(x,0)  +  yF(x,l). 

By  definition,  the  Boolean  sum  of  two  projectors 
(Idempotent  linear  operators)  Is:  P  •  P  * 

P  +  P  -  P  P  .  For  the  time  being, *ve  sKall  assume 
tiiat  tKe  prLsYtlve  function  F  Is  continuous  at  the 
four  comers  of  the  unit  square.  In  which  case  the  pro¬ 
jectors  commute: 


P  F  [FJ 
z  y 


FyP^[Fl . 


(The  main  results  of  this  paper  are,  as  we  shall  soon 
discuss,  concerned  with  problems  Involving  corner  sin¬ 
gularities;  In  those  cases,  the  relevant  projectors  co 
not  cosanute.)  In  terms  of  these  commutative  pro¬ 
jectors  P^  and  Py,  expression  (1.1)  for  U  Is  simply: 

0  •  (Fj  •  Py)[F]  -  (Py  •  P^)(F1.  (1.6) 

Interpolation  schemes  of  this  type  are  known  by 
several  aliases  including  Boolean  am  inter-iolation, 
tranafinite  interpolation,  and  blending  function 
interpolation.  Of  these,  the  term  transflnlte  cones 
closest  to  conveying  the  essence  of  this  class  of 
techniques.  These  methods  are  distinguished  from 
classical  finite  dimensional  Interpolat .on  schemes  by 
the  fact  that  they  Incorporate  a  nondenut-jrahle  -.'inber 
of  scalar  samples  of  F  into  the  interpolant.  More 
precisely,  interpolation  schemes  of  this  class  extract 
from  the  bivariate  primitive  function  F  univariate 
aamplea  of  F,  not  simplv  scalar  samples.  (Note  that 
P  [F]  and  P  [F)  individually  and  •  P  )1F!  are 
transflnite^intcrpolants,  the  product  P  PJF)  • 

P  P  (Fl  Is  merely  the  standard  four  parameter  bilinear 
liiterpolant  to  the  four  corner  values  cf  F.) 

Previous  studies  by  Cordon  and  K.tll  [3j  and  Cordon 
snd  Kelly  [b]  have  been  aimed  at  demonstrating  how,  tnr 
continuous  boundary  conditions,  the  transflnlte,  bi- 
llncarly  blended  interpolant  (1.1)  can  be  employed  to 
reduce  the  computational  effort  in  obtaining  numerical 
solutions  to  second  order  elliptic  boundary  value 
problems.  In  those  two  papers,  the  authors  discuss 
the  following  general  approach:  First,  use  (1.1)  to 
construct  U,  vhlcli  exactly  siaiches  the  given  Dlrlchlet 
boundary  conditions  and  thus  reduces  the  original  pro¬ 
blem  to  one  with  homogrneoue  boundary  condit'ont; 
then,  examine  the  original  problem  for  any  additional 
Information  which  may  bs  inferred  about  the  solution. 
Such  auxiliary  know) edge,  althmigii  perhaps  mcrclv 
qualitative  or  heuristic,  can  often  be  used  to  advan¬ 
tage  In  improving  upon  the  first  ipproxlm.it lou  U 
obtained  by  simply  matching  the  boundary  conditions. 

As  an  example,  the  solution  to  lapl.icc's  equation 
must  satisfy  a  Hixlmum  Principle.  If  the  Initial 
astlsMte  U  does  not,  then  there  arc  simple  ways  ot  con- 
•tnictlng  functions  V  which  vanish  on  3S  snd  arc  such 


that  U  -f  V  both  matches  the  given  boundary  conditions 
and  satisfies  the  Haxinun  principle,  cl.,  [S}-  and  [6]. 
For  the  Poisson  equation  V~U  •  o,  the  sign  of  o  deter¬ 
mines  that  the  solution  is  (locally)  either  sub- 
harmonic  or  superhamonlc  [8],  and  this  auxiliary 
Information  can  be  built  Into  the  exact  boundary 
matching  approximation  U  V.  In  the  actual  testing 
of  these  Ideas,  ve  have  found  Chat  interaotive  com¬ 
puter  graphics  Is  an  almost  Indispensable  aid. 

Numerical  experiments  using  the  techniques 
suggested  In  (5]  and  [6]  to  obtain  good  first  approxi¬ 
mations  with  which  to  enter  standard  iterative  linear 
system  solvers  demonstrated  Chat  very  substantial  re¬ 
ductions  In  total  computational  cost  can  be  realized 
using  such  preprocessing  methods.  Inasmuch  as,  lack¬ 
ing  any  previously  computed  results,  the  standard 
Initialization  of  an  Iterative  scheme  for  solving 
large  linear  systems  is  to  set  all  unknowns  equal  to 
aero  (or  some  constant),  it  is  no  surprise  Chat  an 
exact  boundary  matching  function  which  also  Incorpor¬ 
ates  readily  available  auxiliary  Information  should 
produce  a  more  rapidly  convergent  numerical  solution. 

What  may  be  more  surprising  are  the  results  re¬ 
ported  by  Mitchell,  Marshall  and  Halt  ([9],  [10, 
pp.  174-175J),  and  by  Rice  [11].  Namely,  that  merely 
by  reducing  an  elliptic  problem  with  inhomogeneous 
boundary  conditions  to  a  problem  with  homogeneous  con¬ 
ditions  one  is  able  to  achieve  a  numerical  solution  of 
specified  accuracy  using  far  fewer  algebraic  unknowns. 
(Exact  a  priori  matching  of  boundary  conditions  Is 
tantamount  to  reducing  the  original  problem  to  a  pro¬ 
blem  whose  solution  —  the  "residual"  —  must  satisfy 
homogeneous  boundary  conditions.)  Rice  has  observed 
this  empirically  for  the  collocation  codes  In  the 
ELLPACK  suite,  and  Marshall  and  Mitchell  have  reported 
this  to  be  true  In  experiments  contrasting  standard 
bilinear  finite  elements  with  "exact  boundary, 
elements".  For  the  potential  flow  problem  (7“U  •  0) 
with  a  source  at  (.437,-k),  the  exact  solution  of 
which  is  U  •  log  r  where  r^  •  (x-.437)^  +  (y+k)^, 
Marshall  and  Mitchell  obtained  results  indicating  that 
for  a  "weak"  singularity  at  (.437,-3),  more  than  256 
standard  bilinear  elements  are  required  with  inhomo¬ 
geneous  boundary  conditions  to  achieve  the  same  (four- 
figure)  accuracy  as  can  be  obtained  with  16  elements 
if  the  boundary  conditions  are  first  homogenized.  If 
the  singularity  is  located  at  (.437,-.!),  the  com¬ 
parison  is  roughly  236  elements  to  achieve  three- 
figure  accuracy  with  Inhomogeneous  conditions  versus 
64  with  homogeneous,  cf.  Table  4,  p.  175  of  [10]. 

In  brief,  the  development  of  methodologies  and 
associated  software  preprocessors  to  a  priori  exactly 
match  rather  general  boundary  conditions,  and  thus 
permit  their  homogenization  prior  to  discretization 
and  numerical  solution,  promises  considerable  savings 
in  total  computational  cost,  whether  the  discrete 
linear  system  Is  solved  by  iterative  or  direct 
methods . 

II.  Transflnlte  Interpolation  to  Mixed  (Consistent) 
Boundary  Conditions 

In  this  section,  we  consider  rather  general 
boundary  conditions  of  the  form  oF  +  S(3F/3n)  -  g  on 
the  perimeter  of  the  unit  square.  These  boundary  con¬ 
ditions  are  to  be  thought  of  as  being  associated  with 
some  second  order  elliptic  boundary  value  problem  and, 
without  further  mention,  we  shall  assume  that  they  are 
such  as  to  guarantee  that  the  problem  Is  well-posed. 

In  particular,  this  means  that  the  solution  must  be 
"pinned”  along  at  least  one  of  the  four  edges,  l.c., 
on  at  least  one  of  the  edges  the  function  value  Itself 
must  be  specified. 

In  [6],  Cordon  and  Kelly  considered  mixed  linear 
boundary  cenJitione  of  the  form: 


“  “n^CO.y)  +  B-F  (0,y)  -  e  (y)  along  x-0 


Zii[F]  -  Oj^F(l.y)  •+  “  &i(y)  ’  “long 

MjjlF]  -  IqTCx.O)  +  BqF  (x.O)  -  hgfx)  along  y-0 


(2.1) 


MjlF]  -  o^F(x,l)  +  6^ry(x,l)  -  hj^(x) 


along  y>l 


In  which  the  a  ,  a  ,  6^^  and  are  constants,  and  the 
boundary  conditions  are  consistent,  l.e.: 

ijTfjlF]  -  .VjLj[r]  (l.J  -  0,1).  (2.2) 

In  the  case  of  Dlrichlet  conditions,  the  linear 
operators  and  ere  just: 

£o[F]  -  F(0,y)  -  gp(y),  Lj[F]  -  F(l,y)  -  g^(y) 

•  F(x,0)  -  hQ(x),  Mj^lF]  -  F(x,l)  -  h^(x) 

and  the  consistency  requirement  simply  means  that  the 
boundary  conditions  are  continuous  at  the  four 
comers: 


(2.3) 


WjlBj(y)]  -  ij^[hj(x)]  (1,J  -  0,1).  (2.4) 

Theorem  (Cordon/Kelly)  :  Let  the  and  .V  be  as  In 

(2.1)  and  define  two  projectors  P  and  P  ^as  follows: 

x  y 


FjjIF]  -  6o(x)Ljj[F]  +  6j(x)i^(F] 

P  IF]  -  <^(j(y).’fQ[F]  +  X.j(y)Mj^[F]. 


(2.5) 


where  the  functions  64  end  1I1,  satisfy  the  cardinality 


conditions: 


J 


(Kronecker  Delta)  (2.6) 


for  i,k  -  0,1 

for  j,i  -  0,1. 

Then,  the  function  U  obtained  from  the  Boolean  sum  of 
P  and  P  exactly  satisfies  all  of  the  specified 
boundary^conditions . 

Proof:  The  function  U  is  given  by 

0  -  (Px  +  Py)[F] 

-  ♦p(x)Lp[F]  +  4^(x)Ij[F]  +  <.o(y)Mo[F]  +  :|i^(y).Vj [F] 

-  ♦^(x)<-Q(y)i;^/'p[F]  -  6j^(x)i|.^(y)ij^Wj^[F]. 

The  proof  consists  of  a  straightforward  verification 
ol  the  facts  that  £^[0]  ■  gj(y)  and  .‘■f,  [U]  •  h, (x). 

We  have,  for  example:  ^ 

LoW  -  i^lF)  +  V.o(y)£o.Vo[F)  +  *i(y)Lo.Vj[F] 


(2.8) 


- 

-  8o(y) 

In  which  we  have  used  the  cardinality  conditions  (2.6) 
To  show  that  W,  [U]  “  [F]  •  h  (x)  ,  we  also  use  the 

consistency  hypotheses^ (2. 2) .  O.E.D. 

As  an  Illustration  of  this  result,  consider  a 
function  F  such  that : 


-1.2 


(-1.2)^  +  (y-.9)^ 


ijEFl  -  6P(l,y)  -  F^d.y)  -  6Ilni4-.2)^  +  (y-.9)^+  1.51 

4.  - 

(-.2)^  +  (y-.9)^ 

^^olFl  -  F(x,0)  +  2yy(x,0)  -  ln/(x-1.2)^  +  (-.9)^ 

(2.9) 

1-8  +  1  5 

(x-1.2)-'  +  (-.9)^ 

^llF]  -  F(x.l)  -  la/cx-1.2)^  +  (.1)^  +  1.5 


It  may  be  easily  verified  that  the  four  functions: 


♦o<*>  ■  *  ■  I  ♦i^*^  ■  i 

♦o(y)  -  y  -  1  (ijCy)  -  2  -  y 


(2.10) 


satisfy  the  necessary  cardinality  conditions:  Simply 
apply  the  formulas  for  the  L.  and  M.  to  these  uni¬ 
variate  "blending  functions"  and  coAflrm  the  Kronecker 
delta  properties.  The  corner  values  are  i..V.  [F]  => 
-1.387,  £oM,[Fl  “  -.828,  =  -2.962,  and 

ilWjlF]  “  4.013.  Thus,  the  function  U  given  by 


ll(x,y)  -  (.833-x)- 


1.2 


1.44  +  (y-.9)'‘ 


+  In/.  04  +  (y-.9)^ 


- - 5-  +J  (y-1)  ln/(x-1.2)2  +  . 

3(.04  +  (y-.9)'')  L 

- H -  +  (2-y)ln/(x-1.2)^  +  0.1 

(x-1.2)'  +  .81 J 


+  1. 387 (x-. 833) (y-1)  +  .828(x-.833) (2-y) 

-  1.994 (y-1)  +  .831 (2-y)  +1.5  (2.11) 

satisfies  the  boundary  conditions:  i, [U]  ”  i. [Fj  and 
W,IU)  •  W  (F)  (i,j  =  0,1)  on  the  perimeter  of^S.  The 
graph  of  0  Is  displayed  In  Fig.  2a,  and  Its  contour 
plot  is  depicted  In  Fig.  2b. 

In  the  earlier  vork  by  Cordon  and  Kelly,  the 
authors  assumed  the  existence  of  blending  function 
9,  (x)  and  i|i,(y)  which  satisfy  the  requisite  cardin¬ 
ality  conditions  (2.6).  Here  we  show,  for  boundary 
conditions  of  the  general  form  (2.1),  how  to  actually 
eonstruat  the  i)j(x)  and  ^,(y).  In  particular,  we 
show  that  for  any  choice  bf  the  eight  parameters  o, 

8,,  S,  and  8,  In  (2.1),  we  can  always  find  polynoiuals 
of  degree  three  or  less  which  satisfy  (2.6). 

Lmna:  Let  and  M  ,  and  the  projectors  P  and  P 
be  defined  as  In  the'^ above  theorem.  Then,  ^ere  eilst 
polynomials  of  maximal  degree  three  such  that  (2.6) 
holds. 

Proof:  We  need  carry  out  the  proof  for  only  the  b,(x), 
since  the  F,(y)  are  constructed  Independently  and 
analogously:  To  this  end,  suppose  that  Is  cubic  In 
x:  ” 

♦jj(x)  "  e^j  +  bjjX  +  CjjX^  +  djjX^,  '  (2,12) 

■y  applying  the  linear  operators  and  L.  to  6^ 
and  collecting  terms,  vc  obtain  the  linear  system 

"  "o'o  *  Vo 


-  1 

(2.13) 


’  “l“o  +<V®1>‘*0  +t»i«Bi''-i  +(»i+3Bi)do  -  0 

tor  the  determlnatfon  of  the  polynomial  coefficients 
Sq,  b.,  Cq  and  d^.  (Bear  In  mind  that  the  constants 
and  6^;  whlch'^completely  characterize  and  Lj , 
are  known.)  CleaLiy,  since  there  are  four  unknowns 
and  only  two  equations,  this  system  Is,  In  general, 
underdetermined.  Our  criterion  for  selecting  one 
among  the  (In  general)  two-parameter  family  of 
solutions  Is  to  take  that  solution  which  corresponds 
to  the  minimal  degree  polynomial.  This  resolves  all 
amblqulties  except  two: 

1.  If  aQ(a.+Bj^)  -  “  0  and  both  of  the  following 

2x2  submatrices  are  nonsingular: 

po  °  1  r  °  1 


In  this  case,  we  solve  for  a^  and  Cg  from 
these  and  set  b.  «  dg  “  0. 

2.  If  all  2x2  submatrices  are  singular  es 

h  °  1  f'”  °  1 


•■he  first  of 
:aept 

(2.ii; 


then  solve  for  Sg  and  dg  from  the  first  of  these  and 

•’0  "  *=0  "  °- 

Thus,  assuming  that  cubic  blending  functions  do 
exist,  the  above  procedure  produces  a  unique  bg.  The 
function  ia  obtained  In  a  completely  analogous 
fashion. 

Now,  we  must  show  that  the  linear  system  (2.13) 
td.ll  always  have  at  least  one  solution.  To  see  this, 
consider  the  conditions  under  which  all  five  of  the 
nontrivial  2x2  submatrices  of  (2.13)  are  singular: 


«oC«i+Bj)  -  a^Bg  -  0 

Og(oj+2Bj)  -  0 

OgCoj+BB^)  -  0  (2.16) 

Bg(aj+2Bj)  -  0 

Sg(aj+3Bj)  -  0. 

The  key  to  the  proof  of  existence  is  the  recognition 
that  if  Oi  “  0,  then  B,  cannot  be  zero,  and  vice  versa 
(i"0  and  1).  With  this  in  mind,  it  is  easy  to  show 
chat  the  five  equations  of  (2.16)  cannot  all  hold 
simultaneously.  Q.E.D. 

As  an  example,  consider  the  following  consistent 
boundary  conditions: 


igtFl  -  F(0,y)  +  F^(0,y) 

-  F(l,y) 

Wg[Fl  -  2r(x,0)  +  Fy(x,0) 
M^IF)  -  2F(x,l)  -  Fy(x,l) 


1.5ire"’^  +  ,2 

.25e~^^'*^^cos(ity)  +  .2 

(2.17) 

.25e  (l+sin6vx)  +  .4 
.75e~^^**^(-l+8ln6vx)  +  .4. 


Here,  a  -1,  B.-l.  o  -1,  B  -0,  5-2,  gg-l,  5-2  and 
6,--l,  so  that,  by  Following  the  above  algorithm,  we 
obtain  for  the  blending  functions: 


♦g(x)  -  1  -  X*  ♦j(x)  -  X* 

♦g(y)  -  .5  +  y^  ♦j(y)  -  -y^ 


(2.18) 


which  yield  Che  function: 

U(»,y)  -  1.5»(l-x*)e‘’'  +  .25x^(e ’^^‘‘^^coexy) 

+  .25e"*a+sln6nx)(,5+y’)  (2.19) 

-  ,7Se“^^'‘*^(-l+sln6inc)y^  -1.5x(l-x^)  (.54y'} 

+  4.5iie~^(l-x^)y’  -  .25e"^x^(.54y^) 

-  .75e"Vy^  +  .2. 

By  applying  the  four  linear  operators  and  M,  to 
this  last  expression,  it  can  be  confirmed  that^ 

-  £  [FI  and  W  [U]  =  M  [F]  (i.J  -  0.1),  i.e., 

0  doee  satisfy  the  requisite  boundary  conditions. 

The  graph  and  contour  plot  of  0  are  shorn  in  Figs. 

3a  and  3b. 

Ill .  Trans f ini te  Interpolation  to  Inconsistent 
Mixed  Boundary  Conditions 

As  a  practical  matter,  boundary  conditions  for 
elliptic  problems  are  quite  frequently  not  consistent¬ 
ly  specified.  By  this  we  mean  chat,  although  the 
solution  must  be  smooth  (analytic)  Inside  the  pro¬ 
blem  domain,  it  may  and  often  does  have  singularities 
(discontinuities)  on  the  boundary.  An  elementary 
example  .f  this  is  the  textbook  heat  conduction 
problem  of  du.emining  the  equilibrium  temperature 
distribution  in  a  square  plate,  three  sides  of  which 
are  immersed  in  ice  (0‘C)  and  the  fourth  in  steam 
(100*0 .  (Figs.  Aa  and  Ab  depict  the  graph  of  the 
solution  and  its  contour  plot.) 

In  simple  instances  such  as  this  with  Dirlchlet 
boundary  conditions,  the  analyst  faced  with  solving 
the  boundary  value  problem  will  undoubtedly  be  aware 
of  the  singular  behavior  at  two  of  the  corners  since 
It  is  so  conspicuous.  As  a  rule,  however,  inconsis¬ 
tently  tr-ctified  mixed  boundary  conditions  are  not 
easily  spotted.  For  Instance,  suppose  that  a 
Dirlchlet  condition  F(0,y)  «  8g(y)  Is  specified  along 
the  edge  x~0  of  the  unit  square  and  that  along  Che 
edge  y*0  the  Neumann  condition  F  (r.O)  •  hg(x)  is 
given.  In  order  for  the  solutioS  to  he  smoothly  con¬ 
tinuous  at  the  comer  (0,0),  it  is  necessary  chat 
these  two  conditions  be  consistent,  i.e.: 

11m  ^  F(O.y)  -  lim  F^(x.O) 


The  question  of  consistency  or  Inconsistency  is,  of 
course,  even  more  subtle  for  general  houndary  con¬ 
ditions  of  the  form  (2,1)  above. 

Fortunately,  high  quality  numerical  software  for 
solving  elliptic  problems  is  sufficiently  robust  as  to 
be  able  Co  accommodate  even  grossly  inconsistent 
boundary  conditions.  Provided  chat  a  solution  exists, 
by  taking  a  sufficiently  fine  discretization  (and 
perhaps  employing  some  special  tricks),  Che  applied 
analyst  can  normally  obtain  a  solution  to  whatever 
accuracy  desired.  This,  hotrcver,  is  a  computationally 
expensive  procedure  which  can  be  better  handled  hy 
a  priori  taking  cognizance  of  the  anticipated  singular 
behavior  near  corners.  This  is  the  main  goal  of  this 
section. 

With  the  same  notation  as  Che  previous  section, 
ws  now  consider  botmicrii  operators  L.  and  N.  which  do 
not  ootrmitc:  £,A.'  [F]  y  ,7  ij[Fl  (I,)"'  O,!)";  The  non- 
commutntlvlty  of  the  L,  and  7  Is  wfiat  we  mean  by 
inconaisient  boundary  condiiinna.  Wo  still  find  it 
useful  to  consider  the  projerclon  operators  P^  and  P 


et  (2.5)  in  which  the  dj^fx)  and  *.(y)  are  determined 
In  precisely  the  same  way  as  outlined  in  the  previous 
section.  Now,  however, 

P  P  (F)  y  P  P  [F],  (3.1) 

X  y  y  X 

which  has  Che  important  Implication  that  Che  function 
U  ■  (P  ♦  P  )(F]  will  not  satisfy  the  (Inconsistent) 
boundary  conditions.  (Actually,  (P,.  •  P  )[F]  does 
satisfy  the  conditions  that  £»f(F^  S  P^J^F)]  ■  £.IF) 
(1“0,1),  but  does  not  satlsfv*'che’'other  two  con¬ 
ditions  :  M  [  (P^  •  P  ,)  [  K] )  y  [  F] . ) 

In  order  to  deal  with  corner  inconsistencies,  we 
develop  a  class  of  interpolancs  specifically  designed 
for  the  purpose.  Ac  each  comer  (l,j)  we  construct  a 
special  function  13^^  such  chat 

In  words,  for  fixed  i  and  j,  the  function  U,  vanishes 
under  operation  by  any  of  the  six  Umar  functionals 
£.7  and  M  L,  (kyi,  lyj).  When  operated  on  by  £  7, 
or  B  £  ,  tfie  result  is  £^^.7  [F]  or  .7  £j^[F].  ^  J 

respectively.  ^  ^ 

The  case  of  pure  Dirlchlet  boundary  conditions 

^“o  “  ’  “o  “  “i  '  ■  ®i  ■  ■  ®1  ' 

simplest  to  interpret.  Suppose  the  boundary  con¬ 
ditions  are  such  chat 


F(x.j) 


•  lim  r(x,i)  •  £..7,[F) 
,  JC*i.  ^  ^ 


(l,j-0,l)  (3.3) 


F(i,y)  -  lim  F(l,y)  •  7  £  [ 

J  ^ 


and  £,7,(F)  4  7£  [F];  cf.,  for  example.  Fig.  4a. 
The  function  U.^  which  we  shall  construct  will 
satisfy,  for  tne  case  of  Dirlchlet  conditions: 


11m  U.,(x,J)  -  £.7.IF)  -  lim  F(x.j) 
x*l  ^  ^  x*l  I 

lim  U.,(l,y)  -  7.£JFJ  -  lim  F(l,y) 
T*3  ^  ^  Y*i 


(1,J-0,1)  (3.4) 


and  at  the  three  comers  other  than  (1,3),  U.,  will 
vanish.  ^ 

In  Che  general  case,  suppose  for  the  moment  that 
we  have  the  required  functions  U  (x,y)  which  satisfy 
conditions  (3.2).  Let  W  be  equal^to  Che  sum  of  these 
four  corner  functions; 

W(x.y)  ■  Ugjj(x,y)  +  UQj^(x,y)  +  Uj^g(x,y)  +  Uj^^(x,y), 


Clearly,  W  satisfies  the  eight  conditions: 

-  £j7j[F),  A/j£^[W]  -  7j£^[F)  (i.J-0.1). 


From  this,  we  draw  the  important  conclusion  that,  by 
virtue  of  the  linearity  of  the  operators  £^  and  .7^: 


£j7jlr  -  W]  -  0 
7j£j[r  -  W)  «  0 


for  1,3  "  0,1. 


If  the  solution  to  the  original  intert  olatisn  problem 
is  again  denoted  by  U,  we  want  to  represent  U  as  the 
sins  of  W  and  a  yet  to  be  determined  function  V: 


DU.y)  -  W(x,y)  +  V(x.y). 


(3.8) 


Mow,  since  U  la  Co  satisfy  the  boundary  conditions 
i«[Ol  •  £.[F1  and  W  (UJ  -  «,[FJ  (l.j  -  0.1).  we  have 
ffo«  (3.7)  that:  ^  ^ 


-  0 


(l.J  ■  0,1), 


(3.9) 


which  la  to  say  that  the  function  V  eatiafiea  con- 
aiatant  boundary  aonditiona,  as  defined  In  relation 
(2.2).  Therefore,  we  can  actually  construct  V  using 
Che  techniques  presented  In  the  previous  section. 
Referring  back  to  (3,7),  we  have  that: 


V(x,y)  -  •  Py)lF  -  W) 

-  P,[F  -  W1  +  PytF  -  W], 


(3.10) 


the  last  because  of  (3.7). 

In  summary,  we  first  construct  the  function  U., 
for  each  corner  (l,j).  Then,  we  compute  the  derlvco 
boundary  conditions,  ~  W]  and  W,  (F  -  Vf] ,  and  use 

these  In  expression  (3.10)  for  V.  Tne  function 
U  »  W  +  V  will  then  exactly  satisfy  the  original,  in~ 
ooneiatant  boundary  aonditiona:  LAV]  •  £.  IF)  and 

aePsei_Bar««i^ _ ^  1 


•TfjIUl  -  MjlF]  for  l.j 


0,1. 


Ue  shall  now  without  derivation,  display  the 
functions  U, , .  (For  a  more  complete  treatment,  see 
[7).)  For  every  point  (x,y)  In  S  “  [0,1]  x  [0,1), 
define  the  angles  6  (l,j  -  0,1)  Indicated  in  the 


accompanying  figure: 


tJ 


•qq  ■  arctan(J) 


«10  “  «rctan(^  ~  *) 
1  - 


(3.11) 


*01  "  »tct«n(i-l--)  “  srctan(^  ~"^) . 

He  then  define  the  functions  as  follows: 


+  (i-T(eoo))Vof''" 


Ujj(«.y)  -  ♦o(x)*i(y)lT(ej,^)M^£g[F] 


♦  (l-T(e^j))ijjW^lFll 


♦i(x)*o(y) lT(eio)Wgii  1 


+  (l-T(e^g))igMj[F]] 

«jj(*.y)  "  ♦j(x)ii»j(y)[T(ej^j^)/;^Mj^[F] 


+  (i-T(e^j))Wjr^lF)). 

The  functions  d.(x)  and  i|>,(y)  (i,J  '  0,1)  are  the  same 
as  In  (2.6),  and  the  T(6  ^)  must  satisfy  the  following 
conditions:  ^ 

T(e^j)  -  1  at  Ojj  -  0,  T(e^j)  -  0  at  i  (3.13a) 


*13  ■°*'5e;7"^V 


0  at  e 


se 


2  °  »1J 

13 


V 

13  2 

(3,13b) 


(3.14) 


In  the  case  of  Dlrlchlet  boundary  conditions, 
only  equations  (3.13a)  must  hold,  and  they  are  quite 
simply  satisfied  by  taking: 

26,, 


T(eij) 


(1  - 


4  1.3  -  0.1. 


(3.15) 


For  the  more  general  operators  L,  and  M, ,  the  cubic 
function  ^ 

^®11  2  *®11 

T(6ij)  -  (-^  -  1)  (-7^  +  1)  (1.3  -  0,1)  (3.16) 

satisfies  (3.13a),  (3.13b)  and  (3.14)  as  required. 

Figure  5a  shows  a  perspective  view  and  5b  the 
contour  plot  of  the  Interpolant  to  the  mixed  incon¬ 
sistent  boundary  conditions: 


IglF) 


F(0,y)  -  co8h(^l-y))  +  1 

£j[F)  -  F(l,y)  +  F^(l,y)  -  cosh(Y(l-y))sln(y  y) 
F(x,0) 


Afgirj 


(3.17) 


Fy(x,l) 


"  .5x  +  1 

-  0 


where; 

■  1  -  .5x 

■  1 

and 

-  1 

IfgDgiF)  -  3.509 
LgAfjiFj  -  0 

-  0 


♦o«>> 


♦j(x) 


•  Sx 


♦i(y)  -  y 


D^WglFi  -  2.5 


(3.18) 


WgillF] 


Jf^tjlFJ  -  0. 


(3.19) 


Figure  6a  shows  a  perspective  and  6b  the  contour 
plot  of  the  Interpolant  to  the  mixed  Inconsistent 
boundary  conditions; 

iglFj  -  F(O.y)  -  .25sln(»(4y+.5))  +  .7 

tjlF]  •  F(l,y)  -  2y(l-y)co8(«(2y-.25)  +  .2 


(3.12) 


5 


(3.20) 
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ERRATA 


p.  2,  equation  (2.7),  (P  +  P  ) [F]  should  be  (P  »  P  ) [Fl . 

X  y  X  y 


p.  3,  equation  (2.11), 


(y-1)  ln/(:<-1.2)^  +  .81 - ^ - 

(x-1.2)^  +  .81 


should  be  (y-1) 


In 


n/Cx-1.2)^  +  . 


81  - 


1.8 


(x-1.2)  +  .81 


p.  5,  equation  (3.20),  =  F(l,y)  =  2y(l-y)cos(Tr(2y-.25)  +  .2 


should  be  «  F(l,y)  =  2y(l-y)cos(ii(2y-.25))  +  .2 


